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Wood is a multi-scale material exhibiting a complex viscoplastic response. We study avalanches in 
small wood samples in compression. "Woodquakes" measured by acoustic emission are surprisingly 
similar to earthquakes and crackling noise in rocks and laboratory tests on brittle materials. Both 
the distributions of event energies and of waiting (silent) times follow power-laws. The stress- strain 
response exhibits clear signatures of localization of deformation to "weak spots" or softwood layers, 
as identified using Digital Image Correlation. Even though material structure-dependent localization 
takes place, the avalanche behavior remains scale-free. 


Materials do not deform smoothly in general which 
may be seen on many scales. Avalanches in fracture, 
interface dynamics of domain walls and irreversible plas¬ 
tic deformation exhibit power-law statistical properties 
and are thus "scale-free" P, @|- The intermittency of de¬ 
formation or such non-linear response is a signature for 
collective phenomena Q and is important in practice to 
eventually predict the failure. Fracture avalanches are 
established in many materia ls Ijl , 0|, mostly by acoustic 
emission (AE) technique (5|-[l8[7 and analogies are often 
drawn to "laboratory-scale earthquakes". Apart from the 
presence of power laws in the AE statistics, what other 
universal features can be distinguished is not clear. The¬ 
ory suggests Pi why avalanches with the observed statis¬ 
tical properties exist, but hints about the relevance of the 
detailed physics, including the loading mode, the possible 
difference between 2D and 3D materials, the existence of 
a dominating crack and finally the material properties, 
like disorder and rheology. 

Compressional loading is particularly interesting. 
When the material is crushed, the damaged structure 
and "debris" maintain the ability to carry the load via 
the frictional contacts, and deformation eventually local¬ 
izes into shear bands. Such laboratory experiments may 
also serve as scaled-down models of earthquake faults. 
Recent research concentrates on signatures of criticality 
from models and the role of self-organization of defor¬ 
mation and load [H 0. Experimental studies focus 
on testing the anally of brittle fracture of rocks @-0 
and heterogeneo us |p| -0 or porous media 0|-0| with 
earthquakes [U [2J. Data from small scales presents 
statistical laws like those from earthquake activity. An 
important question across the scales is whether (possi¬ 
ble) correlations in crackling noise and deformation help 
to forecast large events or the final failure 1 , 113 , 0 . 

Many porous materials have a complex structure and 
response to compression. In this Letter, we deal with 
wood, a ubiquitous biological material with a cellular 
structure with hierarchical ordering. The wood cell ar¬ 
rangement depends on the annual growth cycle, leading 
to an alternating series of softwood and hardwood lay¬ 
ers 0, 01 • The mechanical properties vary according 
to the species and origins of the original log and size 
of the sample, and the testing way (along the grains. 


etc.) 0, 0|. Here, we show the presence of avalanches 
in wood compression deformation. Acoustic and optical 
measurements together with strain/stress -curves demon¬ 
strate the existence of avalanches with scale-free proper¬ 
ties, in this "non-brittle" material with complex struc¬ 
ture. Earlier studies of AE in wood fracture exist for 
many applied purposes [26I-0I, but not for the general 
physics. 

Experiments - Small pieces of dry pine (Pinus 
sylvestris) wood are cut from the same plank with some 
variations: 45.6 ± 1.2mm, 10.1 ± 0.6 mm and 9.5 ± 0.7 
mm (length, width, height, averaged over 25 samples). 
The compression, in the vertical direction and perpendic¬ 
ularly to the wood grains, is applied by a piston extension 
on Icm^ area of the sample top (Fig. [TJi, along the grain 
compressions were not tested). An Instron Electropuls 
ElOOO testing machine is used to apply up to 1 kN at 
either constant strain- or stress-rate, of 0.10 ± 0.01 %/s 
and 9 ± 1 kPa/s respectively, the variance in these (slow) 
rates being due to the different sample dimensions. To 
follow the microcrack dynamics, acoustic emission is de¬ 
tected with a piezoelectric transducer, attached on the 
sample upper side. The sampling frequency is 200 kHz. 
After thresholding, events are formed from data, when 
silent times shorter than 35 ps are neglected, in standard 
practice. The event occurrence time t^ is defined as the 
time of the event maximum amplitude and the energy 
Ei is calculated as the square integral of the amplitude 
over the duration. Thus AE time series consist of pairs 
of inter-event waiting time and energy (r^ = U — ti_i, Ei) 
and is analyzed with recent techniques. The number of 
events in a time interval divided by the interval length de¬ 
fines the event rate r{t). Each sample produces typically 
a few thousand events: 4500 ± 2200 and 5200 ± 3300 
with constant strain- and stress-rates, respectively. The 
deformation field is captured by non-invasive digital im¬ 
age correlation (DIG) technique 0, [23,0. Pictures of 
one sample side are taken every second using a grey-scale 
camera (Dalsa Genie HM1024) and are then exploited in 
pairs to obtain the 2-D local strain field 0, . 

Results - Figures Ib-ld demonstrate that AE activity 
or avalanches are present in wood compression. This in¬ 
termittency exists in all stages beyond the linear early re¬ 
sponse: consequent plastic deformation as indicated by a 
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Figure 1. a) Experimental set-up: One end of a wood piece is 
compressed in the machine by the piston extension stick. A 
piezoelectric sensor attached to the other end detects acoustic 
emission. A camera takes pictures of one sample side, enlight¬ 
ened with a LED lamp, b) Compression strain-stress curves 
for two typical experiments at constant strain- and stress-rate 
with AE energies represented with color-code; the curve for 
constant strain-rate test is shifted by +2MPa to avoid curve 
overlap, c) AE events (left y-axis) and cumulative AE energy 
normalized by the total energy (right y-axis) against time for 
constant strain rate compression, d) AE event rate and to¬ 
tal number (left and right y-axis respectively) against time, 
during the same experiment. 


plateau in the stress-strain curve, and stiffening of wood, 
typical of cellular solids, are mixed with avalanche activ¬ 
ity on all scales (Figure [TId), to which question we return 
below for the consequences and origins of the varying 
event rate. For both loading protocols, no large events 
are visible initially, but the activity varies later, as ex¬ 
pected from an avalanching system from small to very 
large events (Figure [TJ: and Id for interesting quanti¬ 
ties: event energy, cumulative AE energy, event rate, 
and the cumulative number of events). Thus come three 
questions, two of them ubiquitous and one particular to 



Figure 2. a) Strain rate, calculated from sample height com¬ 
pression from testing machine and b) AE event rate, both 
for one constant stress-rate experiment; 4 particular instants 
numbered: 3 where high peaks in strain- and AE event- 
ratesare simultaneous, and one reference, c) Strain field from 
Die for the 4 previously defined instants, with fixed color- 
scale for the vertical strain: the localization of deformation 
in softest layers occurs with high AE events- and strain-rates. 
The original image area is also shown (left). 


wood: i) what are the statistical properties of the crack¬ 
ling noise? ii) what are the correlations in the AE signal 
like? iii) does the particular structure of wood become 
important or interact with the scale-free nature of acous¬ 
tic emission or the deformation collective dynamics? 

Localization of deformation and AE - Wood mechanics 
is influenced by the multiscale structure, i.e the annual 
rings. High energy acoustic events are simultaneous with 
rapid reductions in sample height. Acoustic activity is 
also observed outside these periods of high strain rate, 
but the spikes in event rate correspond qualitatively well 
to those in strain rate. Figure [2] shows the mesoscopic 
strain rate et(t) (calculated from the decrease in time of 
the sample height, given by the displacement of the ma¬ 
chine piston) and the event rate r(t) to demonstrate their 
correlation with the localization of deformation from DIG 
[ 2 ^ . The main observed phenomenon is the (partial) col¬ 
lapse of the annual rings ( 2 ^, leading to periods of high 
acoustic activity. Strain maps calculated by the DIG al¬ 
gorithm over particular 10 s time intervals (three over 
intervals of high acoustic activity and one reference) il¬ 
lustrate how the strain localizes at these events. The 
first shows that the strain is localized in two bands cor¬ 
responding to annual rings. The same behavior appears 
in the next two strain maps, but is absent in the last. 
These structure-dependent correlations may be related 
to e.g. the clear precursor events of high energy seen 
right before the major collapses in silica-based porous 
samples (l8| and such rises in event rates are indicative 
of upcoming loss of structural stability Q . 

Statistical properties of acoustic emission - AE data 
from wood compression turns out to follow several em¬ 
pirical deformation or seismological laws. The AE ener- 
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Figure 3. Distribution function of energies of AE events, av¬ 
eraged over all the experiments realized at constant strain 
rate et or stress rate at respectively, exhibit a power law 
P{E) ~ E~^ with P « 1.4, the scaling region and expo¬ 
nent being almost identical in all experiments. Thresholding 
the data with the AE event rate r has little effect on the 
distribution shape. 


gies follows the Gutenberg-Richter (GR) law P{E) 

|31L l82l | for both experiments with a constant strain- and 
stress-rate, with little difference: the P{E) is well fit¬ 
ted, using maximum likelihood method, with a power law 
with exponent 1.37 ± 0.01 for ej constant and 1.43 ± 0.01 
for at constant (Figure [3]). These values are close to the 
ones found in brittle materials such as Vycor (/3 = 1.39 
± 0.02 0 or 1.40 ± 0.05 0 for tests at different stress 
rates), chipboard panels (1.51 ± 0.05 [§|), and for various 
silica-based porous materials ( /3 between 1.35 and 1.55 
0). The exponent dependence on the fitting range is ro¬ 
bust in a large region of around five decades (E from 10“^ 
to 10“^), for both types of experiments, the only change 
being a very slight decrease in high energy event proba¬ 
bility in constant strain rate experiments. In comparison 
to earthquake activity, the rate of acoustic events varies 
considerably over the experiment duration and therefore 
the process is non-stationary. This is taken into account 
by restricting the analysis of AE energy distributions to 
time intervals with event rate above some pre-determined 
threshold value, but these thresholds have almost no ef¬ 
fect on the distributions (figure IS]) . 

Here the waiting time distributions can be roughly di¬ 
vided into two parts depending on the loading protocol 
described by two power laws with different exponents for 
long and short silent times (Figure |4]), arising from large 
variations in the event rate during an experiment (cf. 
the data thresholded with r showing how a cut-off be¬ 
comes visible as expected). At constant strain rate, the 
exponents are -0.97 ± 0.01 and -2.42 ± 0.03 for the two 
regimes divided with waiting times from 0.2 /is to 0.2 s 
and from 0.2 s to 2 s respectively. For constant stress 
rate, the first region is defined for waiting times from 
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Figure 4. Waiting times of AE events are power law dis¬ 
tributed for short waiting times (r ~ 10 “®- 10 “®s) with an 
exponent 7 ~ 1. Considering all data, a second scaling region 
for long waiting times is seen, with 7 « 2. Thresholding by 
the AE event rate r introduces a cut-off to the distributions. 

0.01 s to 1 s with a power law exponent 1.03 ± 0.01 and 
the second exponent is 1.95 ± 0.02, for times between 1 
s and 1000 s. The transition between regions is slightly 
smoother in the experiments with a constant stress rate. 
The difference in the tail exponents (about 2.4 vs. about 
2) should be due to different event rate distributions in 
the two cases; in constant at case, the distribution also 
extends to larger waiting times. The fit with two power- 
laws with exponents around -1 and -2 for short and long 
waiting times respectively is analogous to other porous 
materials 00. The fat waiting time distributions im¬ 
ply temporal correlations in the deformation dynamics 
[ 3 , which we next study by using geophysics techniques 
and recent methods for laboratory experiments. 

Aftershocks and Omori’s law - In the case of earth¬ 
quakes, a mainshock is followed by smaller aftershocks, 
implying the famous seismology Omori’s law [i,0[3i- 
[ 3 ^ : the number of aftershocks per unit time interval 
n{t) is n{t) = K{t + c)“® with K and c constants. The 
modified law relates the aftershock sequence to the time 
elapsed from the mainshock time (l3 | Ims a power law: 
T^Asit—tMs) oc (t—tMs)“^ where r(t) is the event rate of 
aftershocks and p ~ 1. In fracture experiments, events of 
any energy could be considered as mainshocks to verify if 
such statistical laws hold @,0. Thus we vary the main- 
shock energy (threshold), so that all events between two 
subsequent mainshocks are aftershocks of the first one. 
The rates are then averaged over each mainshock energy 
interval. Finally these averaged rates are averaged over 
the different experiments. The event sequences exhibit an 
Omori’s law which becomes less clear as the mainshock 
energy is increased. The Fig. [S^ shows the productivity 
law 0 , 01 , which corresponds to the averaged Omori se¬ 
quences where the aftershock event rate tas is scaled by 
data collapse is obtained for a ~ 0.4. 
The plots for the constant stress rate experiments show a 
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Figure 5. a) Productivity law for different mainshock en¬ 
ergy intervals. The scaled event rate decays with an expo¬ 
nent of 0.75. Empty symbols correspond to constant strain 
rate and filled symbols to constant stress rate.b) Number A of 
AE events for event energy against the waiting time after the 
considered event, in the constant strain rate case; the other 
one is quite similar, c) Number D of events for uncorrelated 
data, obtained by randomizing the (r, i5)-pairs. d) The rela¬ 
tive difference between the original A and randomized D data 
sets shows the presence of correlations, e) Relative difference 
when considering waiting times before an event, with B the 
number of these events. 


a more robust scaling (especially after {t — tMs)= 0.1 s). 
The event rate decays with an exponent around p = 0.75, 
as in Vycor 0 - 

Event energy and waiting time correlations - As shown 
in Fig. [^5, the events are concentrated around certain 
values for interevent times and energies 0: here around 
{TAfter,E) = (0.04 s, 2-10“^). In experiments at constant 
stress rate, this occurs at a slightly lower waiting time, 
but the features are otherwise similar. The waiting times 
after an event of a certain energy present a power law 
relation for aftershocks (high energies followed by short 
waiting times) 0[3ii, which spans around five decades 


in energy with an exponent around 4. This relation is 
absent in the waiting time before an event. Some cor¬ 
relation between small energies and small waiting times 
is present instead. The comparison of the experiments 
at constant strain rate with an artificial uncorrelated 
dataset (Fig. (S]:) (the (r, £l)-pairs randomized, and re¬ 
peated 100 times for better statistics) shows a gap with 
less events with energies in the most prominent range 
from 10“^ to 10“® with waiting times from around 10“^ 
to 10“^ s. A similar but much smaller gap is found for 
constant stress rate. The relative differences between 
waiting times after an event and uncorrelated waiting 
times show that aftershock correlations exist (Fig. [5ji), 
more clearly than in porous material like Vycor [13| ■ The 
relative differences between the foreshock waiting times 
and uncorrelated waiting times show a correlation be¬ 
tween small energies and small waiting times (Fig. [5^), 
and a slight correlation between high energy and high 
waiting time events across all energies and waiting times 
over 10“^ s. 

Conclusion - Wood compression shows clear evidence 
of scale-free avalanche activity, both in stress- or strain- 
rate controlled protocols. The statistics of the avalanches 
(distributions of AE events energies and waiting times) 
exhibit exponents reminiscent of brittle porous materials 
and geophysics. The dynamics of compression is corre¬ 
lated in further analogy with rocks and other such materi¬ 
als and indeed even earthquakes, as demonstrated by the 
modified Omori’s law, the law of aftershock productivity 
and the universal scaling law for the waiting time dis¬ 
tribution typically used in statistical seismology. Scale- 
free avalanches co-exist in wood in spite of the material 
structural features, clearly visible and relevant in com¬ 
pression curves or by microscopic analysis of the larger 
events (DIG). We conjecture that this fact is only man¬ 
ifested in terms of a varying event (or production) rate. 
Further studies would be needed to analyze whether soft¬ 
wood layer collapse is always correlated with an increased 
AE activity, and whether such large events are in fact 
scale-free or imply temporary changes in acoustic activ¬ 
ity correlations. Eurther work is needed to advance mod¬ 
els to reproduce these results, including also the wood 
cellular structure. Similarly, experiments are called for 
to explore the role of rheology (temperature [l^, [H, 0 
and humidity [IJ, 0|) in wood and in other candidate 
materials of natural origin. 
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